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Abstract It is well known that many real world networks have the power-law degree distribution 
(scale- free property). However there are no rigorous results for continuous-time quantum walks on 
such realistic graphs. In this paper, we analyze space-time behaviors of continuous-time quantum 
walks and random walks on the threshold network model which is a reasonable candidate model 
having scale-free property. We show that the quantum walker exhibits localization at the starting 
point, although the random walker tends to spread uniformly. 

1 Introduction 

Continuous-time quantum walks, which are the quantum counterparts of the classical ran- 
dom walks, have been widely studied on various deterministic graphs, such as the line |17] . 
star graph [33]H0] . cycle graph [2|I15P29] . dendrimers [28], spidernet graphs [35], the Dual Sier- 
pinski Gasket [T], direct product of Cayley graphs [SB], quotient graphs [52], odd graphs [55] . 
trees [TB|[TS] and ultrametric spaces [IH]. For further information, see reviews such as [^DIES]- 
Also there are simulation based study of continuous-time quantum walks on probabilistic 
graphs, such as small- world networks and Erdos-Renyi random graph [41j. However there 
are no rigorous results for continuous-time quantum walks on such probabilistic graphs. In 
this paper, we focus on the continuous-time quantum walk on a random graph called the 
threshold network model. 

Many real world networks (graphs) are characterized by small diameters, high cluster- 
ing, and power-law (scale-free) degree distributions [3llll|3T]. The threshold network model 
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belongs to the so-called hidden variable models [71|38] and is known for being capable of 
generating scale-free networks. Their mean behavior [SlITlIinillllESlEZlES] and limit theo- 
rems P[I2l[I3l[2I] for the degree, the clustering coefficients, the number of subgraphs, and 
the average distance have been analyzed. The strong law of large numbers and central limit 
theorem for the rank of the adjacency matrix of the model with self-loops are given by [6]. 
Eigenvalues and eigenvectors of the adjacency matrix [14J, the Laplacian matrix [26l[27] of 
the model have been studied. See also [SllI2l[I3l[2Il[22l[24, 25j for related works. 

This paper is organized as follows. We define the threshold network model and give a 
brief review of the hierarchical structure of the graph in Section 2. In Section 3, we define 
the continuous-time quantum walk on the threshold network model and the special setting 
called the binary threshold model. The main results are presented in this section and the 
proofs are given in Section 4. Results on the continuous-time random walks on the models 
are obtained in Section 5. Summary is given in the last section. 

2 Threshold network model 

The threshold network model Qn{X, 6) is a random graph on the vertex set V" = {1, 2, . . . , ri}. 
Let {Xi, X2, . . . , Xn} be independent copies of a random variable X with distribution P. We 
draw an edge between two distinct vertices i,j G V if Xi + Xj > 6 where 6' G M is a constant 
called a threshold. Hereafter, we use P°° as the distribution of {Xj}^^. 

Each sample graph G G Qn{X,d) has a hierarchical structure described by the so-called 
creation sequence [HIITT] . Let X{i) < X{2) < ■ ■ ■ < X{n) be a rearranged sequence of random 
variables Xi, X2, . . . , X„ in increasing order. If X(i) + X(„) > 6*, we have 



which means that the vertex corresponding to is connected with the n — 1 other vertices. 
Otherwise, we have 



which means that the vertex corresponding to X{i) is isolatedD We set Sn = 1 or = 
according as the former case or the latter occurs. Then, according to the case we remove the 
random variable X(„) or X(i), we continue similar procedure to define . . . , S2. Finally, 
we set Si = S2 and obtain a {0, l}-sequence {si, S2, . . . , which is called the creation 
sequence of G and is denoted by Sq- 

Given a creation sequence Sg let ki and k denote the number of consecutive bits of I's 
and O's, respectively, as follows: 



It may happen that fci = or = 0, but we have /c2, • • • , ^m, ^1, • • • , lm~i > 1 and m > 1. 
Moreover, by definition we have two cases: (a) ki = (equivalently Si = 0) and li > 2; (b) 
ki >2 (equivalently, Si = 1). 



For example, if Sg = {1,1,0,0,1,0,1,0} then ki = 2, Zi = 2, /c2 = 1, /2 = 1, k^ = 
1, /s = 1 and Fig. [1] shows the shape of G. 



9 < X(i) + X(„) < X(2) + X(„) < ■ ■ ■ < X(„_i) + X(„), 



9 > -^(1) + X(^n) > ■ ■ ■ > X(i) + X(3) > X(i) + X(2), 





(2.1) 
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Figure 1: A threshold graph G corresponding to Sq = {1, 1, 0, 0, 1, 0, 1, 0} 



The creation sequence Sq gives rise to a partition of the vertex set: 

m m 
i=l i=l 

The subgraph induced by v}^^ is the complete graph of ki vertices, and that induced by V^^^ 
is the null graph of h vertices. Moreover, every vertex in v}^^ (resp. v}^^) is connected (resp. 
disconnected) with all vertices in 

vl'^u---uv}'^uvl'^u---uv;^\. 

In general, a graph possessing the above hierarchical structure is called a threshold graph [22] . 



3 Our model and results 

Let Ag be the adjacency matrix and Dq be the diagonal matrix of degrees (the sum of the 
rows of Ag) of G E Qn{X, 0). Then the Laplacian matrix Lg of G is given by Lg = Dg — Ag- 
The time evolution operator of a continuous-time quantum walk on G is defined by 

= ^""^ - E -^^o- (3.2) 

A:=0 

Let {'^nt}t>o be the probability amplitude of the quantum walk, i.e., \E'^( = U^^'^^q, 
and Xn^t denotes the position of the quantum walker at time t. Then the probability that 
the quantum walker on G is in position x G ^ at time t with initial condition ^,^q is defined 
by 

P^iXn, = x) ^ \^l,{x)\\ 

where ^^^^ = -^[^'^^^(l) • • • ^^^(n)] . Here '^A denotes the transpose of a matrix A. 
The time evolution operator is obtained as follows: 
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Theorem 3.1 Suppose G G ^„,(X, 6') is connected. The time evolution operator U^^ of the 
continuous -time quantum walk on G is given by 



^fv G V}'\ w G Vj^ U V}'^ {J<^- 1), 
Here A^^^ denotes the {v,w) element of a matrix A and 



{U^tf''^ = ( U^^iw) - \- -) e^*(^^-+i) + V 

\ % I / j=i+ 



m— 1 / T-, n \ HDi . ^ 



Where D^^ and Di^ denote the degree of vertices in V^^"* and V^'^\ respectively and Ia{x) is 
the indicator function of a set A, i.e., Ia{x) = 1 if x & A and Ia{x) = otherwise. 

Theorem 13.11 shows that we can obtain the probabihty of the quantum walker in position 
X E V ai time t for any initial conditions at least in principle. But in general cases, it is hard 
to obtain the probability. In this paper, we analyze behaviors of quantum walks starting 
from a vertex f , i.e., the case of 



1, if s = V, 
0, otherwise. 



Theorem 3.2 Suppose G is connected. The limit of the probability of the quantum walker 
starting from a vertex v G Vm^ is given by 

{1 iij' X —— v 
' ' for -almost every G . 

0, otherwise, 

Because P^^ is not converge in t — )■ cxd, we study the time-averaged probability P^ix) 
defined by 

T 



P^{x) = hm i / P„^,(X„, = x)dt. 
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Theorem 3.3 Suppose G G ^„(X, is connected. The time- averaged probability P^{x) of 
the quantum walker on G starting from a vertex v G Vm^ is 



otherwise. 



In order to study more detailed properties of the quantum walk on the model, we focus 
on the threshold network model Qn{X, 9) defined by Bernoulli trials with success probability 
p e (0, 1), i.e., P(X = !) = !- P(X = 0) = p, and a threshold 6 G [0, 1). We call this 
model the binary threshold model Qnip)- For sach G G Qnip)^ realization G of Gn{p), we 
consider a partition of the vertex set V: 



V 



(1) 



{t:X, = 1}, Vl'^ = {t:X, = 0}. 



G 



It is easy to see that the subgraph induced by Vq^ is the complete graph on kc = ^Vq^ 
vertices, and that induced by Vq^ is the null graph on Iq = ^Vq^ vertices, where is the 
number of elements in a set A. Moreover, every vertex in Vq^ is connected with all vertices 
in V^""*. Note that this is the case of m = 2, ki = I2 = 0, /c2 = and h = Iq in Eq. (12. ip . 

The time evolution operator of the continuous-time quantum walk on the binary threshold 
model is obtained as follows: 

Theorem 3.4 The time evolution operator U^^ of the continuous-time quantum walk on 
G G Gn{p) is given by 



UkG,kG 



UkG,lG 

Uig,Ig 



The elements of U^^ are 



UkG,kG 
UkG,lG 
UlG,kG 

Uig,Ig 



e'^^'lkG + 



1-e 



int 



n 



'-kG,kGi 



n 



Ant 



n 



^kG,lc 



'-lG,kG^ 



e^kGtj^^ + ( - + 



kce 



int 



^ikGt 



n 



nl 



G 



G 



^Ig,Ig^ 



where Ij j is the i x j matrix consisting of only 1 and li is the i x i identity matrix. 

On the binary threshold model, a strong localization is observed at any starting point as 
follows: 

Proposition 3.5 The limit of the probability of the quantum walker starting from a vertex 
V G V is given by 



lim P^^,{x) 



, ifx = v, 
0, otherwise, 



for -almost every G. 
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Note that if the quantum walker starts from a vertex v G V^^^ then the statement of Propo- 
sition 13.51 is the same as Theorem 13.21 

The time averaged probabihty of the quantum walker on the binary threshold model is 
obtained as follows: 

Proposition 3.6 The time-averaged probability Pn{x) of the quantum walker onG E Qn{p) 
starting from a vertex v G Vq^ is given by 



nlr 



4 + '^ 

2_ 

2 ) 



+ 



7,2; 



^/xGri°^M, 
otherwise. 

(1) 



Note that when the quantum walker starts from a vertex f G , we can use Theorem 13.3 



4 Proofs 

The Laplacian matrix Lq of G G Qn{p) is given by 



J^G — UT ■ 

Eigenvalues and eigenvectors of Lg are known as follows pT | [26 | l27]: 



eigenvalue 


eigenvectors 


n 




-J 

_On-j_l,l_ 




kc 




.0'G-i-i,iJ 





W/g = . 



where Ojj- is the « x j zero matrix. Note that the set of these eigenvectors forms an or- 
thonormal basis of M". Thus we can define a orthogonal matrix Bg corresponding to the 
eigenvectors as follows: 



Eg = [vi ... Vfcg wi ... Wi^] . (4.3) 
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Because U^^ is an n x n (finite) matrix, using Eqs. (13 ■2p and (14. 3p . it is represented by 



Un,t = Bg 



Oi 



ka 



Oi 



la-l 



'Bo. 



It is easy to see that 



G 
n,t 



ka nk(. 



Ig \ c'int _|_ 1 



^kckc 



i(l -e" 



int\ 



^ka. 



int\ 



L«G,«G 



By using a relation — l/Zc^ + Ic/nkc = — 1/"^, we obtain Theorem [X 
It is given by simple calculations that 



1 - - 1 e^"* + - 
n / n 



- (1 - e*"*) 



n 



1 1 (1 - cosnt) 



n 



n 



— (2 — cos nt) 



oikat 



+ 



int 



uIq ~'~ n 



1 



G 



+ 



2t 



G 



nl 



G 



1 



^ikQt 



+ 



kae 



int 



l-G 



uIq ~'~ n 



1 

72" 



G 



COS let + -, — COS kat ] + 



G 



2k 



G 



nH 



cos nt, 



G 



1 

+ — 



cos lat + 7^ cos — — cos rat 



On the other hand, we see that lim„_>oo ^g/""- = P and lim„_>oo ^g/""- = 1 — P°°-almost 
surely by the strong law of large numbers for the i.i.d. Bernoulli sequence. Combining these 
facts and Theorem 13 ■4[ we have Proposition 13.51 Also we obtain Proposition [3]6] immediately 
from 



lim - , 



COS nt dt = lim — 



COS kat dt = lim — 



cos let dt = 0. 



In the case of G G Qn{X,6) which is connected, eigenvalues and eigenvectors of Lg are 
also known as follows pT | l26 | l27]: 
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eigenvalue 


eigenvectors 


{2<i<m) 


1 

1 

y/ (ki+di)k 


( 


-J 

di^ki,\ 


(l<J</c^-l), 


Dk 

(2 < i < m - 1) 


1 

1 

^/(ki+di)l, 


( 


1 

-j 

(/Ci + d 

—li^k 


(1< 

fit 

i+di,l 


J<h-l): 




1 


-i 






1 


-J 

.0d2-j-l,l. 


(j = l,...,/l-l), 


Du 

11 


1 

(fci+(ii);id2 


(/ci + di)li^^i 

-ll^ki+di,! 


, ifA;i^O, 
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where Ui = Ylij>A^3 + h) ^^'^ ~ Ylij<A^j + h)- By the same argument in the proof of 
Theorem 13.41 and the following relations: 

h = Dk, - + 1 (if h ^ 0), 

fc, = - A, (2 <2 <m), 

h = Dk^^^-Dk^ (l<2<m-l), 

we have Theorem 13.11 

Note that = n — 1 and Di^ = by assumpution. Comparing Theorem 13.11 with 
Theorem 13. 4[ the transition probability of the walker on G G Qnip) starting from a vertex 

V e Vm^ is the same as that of on G G Qn{j>) starting from a vertex v G V^"'^^ Thus we have 
Theorems 13.21 and 13.31 

5 Classical Case 

The time evolution operator lA^^ of a continuous-time random walk on G G Qn{X,&) is 
defined by 

7/G _ -the, - (~^)^ Tk 

k=0 

Let {Vnt}t>o be the probability distribution of the random walk, i.e., V^t = ^nt^no^ ^'^d 
Yn^t denotes the position of the random walker at time t. By the same observation in the 
previous section, we have the same results for U^^ as Theorems 13.11 and 13.41 by exchanging 
it of t/^j for t. Using these results, we have the following: 

Proposition 5.1 The limit of the probability that the random walker starting from a vertex 

V G Vm^ is given by 

lim nV^tiy) = I5 for all y & V , for -almost every G. 

n—>OD ' 

Proposition 5.2 The long-time limit of the probability of the random walk on G E Qn{p) 
starting from a vertex v E V is given by 

limP5(y) = -, for ally eV. 

We can also estimate the time-averaged probability Pniv)- By simple calculation, we have 
Vniv) = 1/""- ^ random walk on G G Qn{p) starting from a vertex v eV. 

6 Summary 

In this paper, we study the continuous-time quantum and random walks on the threshold 
network model. By comparing Theorem 13.21 with Proposition 15.11 we have quite different 
limit behaviors in n — t- 00 for the two types of walks starting from a vertex which degree 
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equals n — 1. Although quantum walkers exhibit strong localization at the starting point, 
random walkers tend to spread uniformly. 

Theorem 13.31 and Proposition 13.61 show that the time-averaged probabilities of quantum 
walkers are not the uniform distribution (different from random walks). Furthermore, the 
time-averaged probability shows localization at starting point as n — >■ cx3. In the case of 
the binary threshold model, the rate of convergence are slightly different in the two starting 
points. Indeed, we obtain lim„_>.oo n{l — P^{vi)) = 2 < lim^^oo '^(l — Pni'^o)) = 2/(1 — p), 
P°°-almost surely for vq G Vq^ and Vi G Vq\ A study covering the more general setting is 
now in progress. 
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Society for the Promotion of Science (Grant No. 21540118). 

References 

[1] Agiiari, E., Blumen, A., and Miilken, O. (2008) Dynamics of continuous-time quantum walks in re- 
stricted geometries. J. Phys. A 41, 445301. 

[2] Ahmadi, A., Belk, R., Tamon, C, and Wendler, C. (2003) On mixing in continuous-time quantum 
walks on some circulant graphs. Quantum Inf. Comput. 3, 611-618. 

[3] Albert, R., and Barabasi, A. -L. (2002) Statistical mechanics of complex networks. Rev. Mod. Phys. 
74, 47-97. 

[4] Boccaletti, S., Latora, V., Moreno, Y., Chavez, M., and Hwang, D. -U. (2006) Complex networks: 
structure and dynamics. Phys. Rep. 424, 175-308. 

[5] Boguha, M., and Pastor-Satorras, R. (2003) Class of correlated random networks with hidden variables. 
Phys. Rev. E Q8, 036112. 

[6] Bose, A., and Sen, A. (2007) On asymptotic properties of the rank of a special random adjacency 
matrix. Elect. Comm. in Probab. 12, 200-205. 

[7] Caldarelli, C, Capocci, A., De Los Rios, P., and Mufioz, M. A. (2002) Scale-free networks from varying 
vertex intrinsic fitness. Phys. Rev. Lett. 89, 258702. 

[8] Diaconis, P., Holmes, S., and Janson, S. (2009) Threshold graph limits and random threshold graphs. 
Internet Mathematics 5, no. 3, 267-318. 

[9] Fujihara, A., Ide, Y., Konno, N., Masuda, N., Miwa, H., and Uchida, M. (2009) Limit theorems for 
the average distance and the degree distribution of the threshold network model. Interdisciplinary 
Information Sciences 15, no. 3, 361-366. 

[10] Fujihara, A., Uchida, M., and Miwa, H. (2010) Universal power laws in threshold network model: 
theoretical analysis based on extreme value theory. Physica A 389, 1124-1130. 

[11] Hagberg, A., Schult, D. A., and Swart, P. J. (2006) Designing threshold networks with given structural 
and dynamical properties. Phys. Rev. E 74, 056116. 

[12] Ide, Y., Konno, N., and Masuda, N. (2007) Limit theorems for some statistics of a generalized threshold 
network model. RIMS Kokyuroku, no. 1551, Theory of Biomathematics and its Applications III, 81-86. 



10 



[13] Ide, Y., Konno, N., and Masuda, N. (2009) Statistical properties of a generalized threshold network 
modcL to appear in Methodol. Comput. Appl. Probab. 

[14] Ide, Y., Konno, N., and Obata, N. (2010) Spectral properties of the threshold network model, preprint. 
larXiv:1001.ni86l 

[15] Inui, N., Kasahara, K., Konishi, Y., and Konno, N. (2005) Evolution of continuous-time quantum 
random walks on cycles. Fluctuation and Noise Letters 5, L73-L83. 

[16] Jafarizadeh, M. A., and Salimi, S. (2007) Investigation of continuous- time quantum walk via spectral 
distribution associated with adjacency matrix. Ann. Phys. 322, 1005-1033. 

[17] Konno, N. (2005) Limit theorem for continuous-time quantum walk on the line. Phys. Rev. E 72. 
026113. 

[18] Konno, N. (2006) Continuous-time quantum walks on trees in quantum probability theory. Inf. Dim. 
Anal. Quantum Probab. Rel. Topics 9, 287-297. 

[19] Konno, N. (2006) Continuous-time quantum walks on ultrametric spaces. Int. J. Quantum Inf. 4, 
1023-1035. 

[20] Konno, N. (2008) Quantum walks, U. Franz and M. Schurmann (Eds): Quantum Potential Theory, 
Lecture Notes in Mathematics, Vol. 1954, Springer, pp. 309-452. 

[21] Konno, N., Masuda, N., Roy, R., and Sarkar, A. (2005) Rigorous results on the threshold network 
model. J. Phys. A: Math. Gen. 38, 6277-6291. 

[22] Mahadev, N.V.R., and Peled, U.N. (1995) Threshold Graphs and Related Topics, Elsevier. 

[23] Masuda, N., Miwa, H., and Konno, N. (2004) Analysis of scale- free networks based on a threshold 
graph with intrinsic vertex weights. Phys. Rev. E 70, 036124. 

[24] Masuda, N., Miwa, H., and Konno, N. (2005) Geographical threshold graphs with small- world and 
scale-free properties. Phys. Rev. E 71, 036108. 

[25] Masuda, N., and Konno, N. (2006) VIP-club phenomenon: Emergence of elites and masterminds in 
social networks. Social Networks 28, 297-309. 

[26] Merris, R. (1994) Degree maximal graphs are Laplacian integral. Linear Algebr. Appl. 199, 381-389. 

[27] Merris, R. (1998) Laplacian graph eigenvectors. Linear Algebr. Appl. 278, 221-236. 

[28] Miilken, O., Bierbaum, V., and Blumen, A. (2006) Coherent exciton transport in dendrimers and 
continuous-time quantum walks. J. Ghem. Phys. 124, 124905. 

[29] Miilken, O., and Blumen, A. (2006) Continuous time quantum walks in phase space. Phys. Rev. A 73, 
012105. 

[30] Miilken, O., Pernice, V., and Blumen, A. (2007) Quantum tranceport on small-world networks: A 
continuous-time quantum walk approach. Phys. Rev. E 76, 051125. 

[31] Newman, M. E. J. (2003) The structure and function of complex networks. SIAM Rev. 45, 167-256. 

[32] Salimi, S. (2008) Study of continuous-time quantum walks on quotient graphs via quantum probability 
theory. Int. J. Quantum Inf. 6, 945-957. 

[33] Salimi, S. (2008) Quantum central limit theorem for continuous-time quantum walks on odd graphs in 
quantum probability theory. Int. J. Theor. Phys. 47, 3298-3309. 

11 



[34] Salimi, S. (2009) Continuous-time quantum walks on star graphs. Annals of Physics 324, 1185-1193. 

[35] Salimi, S. (2010) Continuous-time quantum walks on semi-regular spidernet graphs via quantum prob- 
ability theory. Quantum Inf. Process. 9, 75-91. 

[36] Salimi, S., and Jafarizadeh, M. (2009) Continuous-time classical and quantum random walk on direct 
product of Cayley graphs. Commun. Theor. Phys. 51, 1003-1009. 

[37] Servedio, V. D. P., Caldarelli, G., and Butta, P. (2004) Vertex intrinsic fitness: How to produce 
arbitrary scale-free networks. Phys. Rev. E 70, 056126. 

[38] Soderberg, B. (2002) General formalism for inhomogeneous random graphs. Phys. Rev. E 66, 066121. 

[39] Venegas-Andraca, S. E. (2008) Quantum Walks for Computer Scientists, Morgan and Claypool. 

[40] Xu, X. P. (2009) Exact analytical results for quantum walks on star graph. J. Phys. A: Math. Theor. 
42, 115205. 

[41] Xu, X. P., and Liu, F. (2008) Continuous-time quantum walks on Erdos-Renyi networks. Phys. Lett. 
A 372, 6727. 



12 



